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ABSTRACT 


a digital feasio er representation Ot the ASQ-81 
magnetometer is derived from the s-plane transfer functions 
of the system through the use of a bilinear transformation. 
A FORTRAN computer program is written which applies’ this 
representation to time-sampled total magnetic field data in 
order to obtain a time series representation of ASQ-81 
exltered total field. A series of simulations and a field 
experiment are conducted which verify the program output. 
Mmaplications of this program include usage in conjunction 
with geomagnetic field data in order to produce a new data 
set representative of geomagnetic noise observed by Navy MAD 
(Magnetic Anomaly Detection) aircraft with the potential to 
investigate techniques of reducing geomagnetic noise in MAD 


eaercraft. 
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I. INTRODUCTION 


The detection and location of submarines (and other 
Magnetic bodies) through the discrimination of changes. or 
anomalies in the Earth's magnetic field is called Magnetic 
Anomaly Detection or MAD. In this technique, a magnetometer 
measures the magnitude of the Earth's magnetic field and 
provides an indication of that magnitude, or, more usually, 
an indication of changes in the magnitude of the Earth's 
field. These changes, or anomalies, can indicate the 
presence of magnetized bodies which may or may not be a 
submarine. 

The magnetometer currently in use in the United States 
Navy for use in this MAD process is the AN/ASQ-81 
metastable helium vapor total field magnetometer. 

Research is currently being conducted at the Naval 
Postgraduate School in Monterey, California, in various 
aspects of the applications of magnetometers, including 
Magnetic Anomaly Detection (MAD). Within the context of 
this research, magnetic field measurements are made through 
the use of sets of wire wound coils vice any specific 
magnetometer or magnetic detecting system. The data 


collected through the use of these coils is evaluated and 





processed in a variety of methods for different project 
goals. 

This thesis project is designed to produce an acceptable 
alternative to the physical presence of an experimental 
AN/ASQ-81 magnetometer at the postgraduate school by 
allowing the determination, in conjunction with other 
research in progress, of the output of the AN/ASQ-81 
magnetometer from the data collected from the _ school's 
measurement coils. It is hoped that this will assist future 
research projects as, for example, in allowing a 
determination of environmental noise of such characteristics 
as to affect the AN/ASQ-81 magnetometer operationally with 
the eventual goal of providing an environmental noise index 
or a system of removing such noise from the magnetometer- 


detection system. 


10 





IIT. GEOMAGNETICS REVIEW 


A. EARTH'S MAGNETIC FIELD 


1. Constituents of the Geomagnetic Field 








The most common method of specifying the consituting 


Memes of the geomagnetic field is to divide the field in 


meme Of distance from the center of the Earth. This 
method results whe! three classifications: intemal, 


crustal, and external. [Ref. 1] 

The internal field originates in the core region and 
is the most stable field, containing only extremely low 
frequency temporal variations. The crustal, or anomalous, 
field arises from modifications made on the internal field 
by materials and structures in the Earth's crust. These 
Variations are not constant with regard EO spatial 
locations, and comprise part of what is known as geological 
Variations. The external field is the most dynamic and 
arises from many sources, including the interaction between 
the solar wind and the Earth's magnetic field. 

In addition to this method of defining the Earth's 
Magnetic field is the method of time variations. This 
method consists of considering that part of the field which 


varies with periodicities greater than about one year as the 








steady field and everything else as the variation field. 
(Ref. 2] 

The steady field consists of the internal field, also 
refered to as the main field. Slow variations of the main 
field with periods of years or longer are referred to as 
secular variations. 

There are various elements that contribute to. the 
geomagnetic field, some of which are external to the Earth's 
surface. External contributions make up only a small part 
of the steady field, but play an important role in the 
variation field. These external sources include current 
systems in the Earth's upper atmosphere affected by solar 
electromagnetic radiation and gravitation, solar corpuscular 
radiation and the interaction of solar plasma with the main 
field, and the effect of the solar inteplanetary fieid. 
(Ref. 3] 

The geomagnetic field changes with time. As 
previously mentioned, very slow variations in the main field 
With periods of on the order of years to thousands of years 
are referred to as secular variations. Secular variations 
are caused by a variation in the strength or orientation of 
the Earth's center dipole. 

Other time variations of the field can be categorized 
into quiet variation fields and disturbed variation fields. 
Disturbed Variation fields include eeomagnetic 


micropulsations, which are of particular interest to 


1 





operational forces as these can mask target signatures’ and 
are therefore a source of noise to MAD sensors. 

Quiet variation fields are those which are not due to 
disturbances in the interplanetary environment and which 
vary slowly and regularly. [Ref. 3] 

Disturbed variation fields are geomagnetic field 


variations that appear to be the result of interplanetary 


environmental changes and do not posess a simple 
Meritodicity. These Variations include ionospheric 
disturbances, the aurora, geomagnetic storms, and 


geomagnetic micropulsations. 
@ Elements of the Magnetic Field Vector 

The geomagnetic field vector is characterized at any 
Dommt Dy its direction and magnitude. This is commonly 
accomplished through a system of coordinates as shown in 
meeure 2.1. The field is measured in terms of local 
coordinates with respect to true North. [Ref. 3] 

The various coordinates are referred to as magnetic 
elements and are defined as follows: 

B: Total field intensity (the symbol F is 


sometimes also used, as in this figure.) 


H: Horizontal component 

X: Northward, or NorthSouth component 
Y: Eastward, or EastWest component 

Z: Downward, or Vertical component 


ie 
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Figure 2.1 Magnetic Field Elements [Ref. 4]. 


D: Declination or magnetic variation 


This is the angle between X and H and 


is measured positive eastward. 


ie inetinataon or Dip Angle. 


This is the angle between H and B (or F) 


and is measured positive downward. 
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III. THE AN/ASQ-81 MAGNETOMETER 


Ae DESCRIPTION OF SYSTEM OPERATION 

The Magnetic Anomaly Detecting set currently in use in 
the US Navy is the AN/ASQ-81 magnetometer. This set is 
used to locate and classify submerged submarines by sensing 
disturbances in the Earth's magnetic field (anomalies) 
caused by the presence of the magnetic mass of the 
submarine. The disturbance of the Earth's field is detected 
by the magnetometer, processed through filtering circuits, 
and amplified. The output signal of the magnetometer is 
displayed on a chart recorder for interpretation by an 
operator. 

The magnetic detecting set is a metastable helium vapor 
magnetometer. The operation of the magnetometer is based on 
the light absorbtion properties of helium gas subjected to 
certain light stimulus (optical pumping), radio frequency 
excitation, and the Earth's magnetic field. The 
Magnetometer consists of a helium lamp, lens and polarizer 
to generate a beam of polarized light radiation. This 
focused and polarized light beam is directed through a 
helium absorbtion cell to an infrared (IR) detector. Some 
of the helium gas in the absorbtion cell is maintained in a 


metastable state by application of VHF excitation. 


19s, 
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The Earth's magnetic field imposes a magnetic force upon 


the excited helium vapor atoms to force the atoms into one 


of three energy sublevels. This is called the Zeeman 
effect. The rate or frequency of atomic precession caused 
by this effect is called the Larmor frequency. A helium 


lamp is used to optically pump the atoms in the absorbtion 
cell, with the result that the polarized light energy 
passing through the absorbtion cell will polarize 
(magnetize) the helium atoms in the absorbtion cells by 
selectively pumping the Zeeman levels of the energy of the 
helium atoms in the cell. The magnetization direction is 
determined by tthe polarization of the photons from the 


helium lamp. 
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RF energy is then introduced to the absorbtion cell in 
the form of an additional magnetic field imposed through the 
use of coils oriented perpendicular to. the precessed 


polarized helium atoms in the absorbtion cell and energized 


by a variable frequency RF oscillator. The RF oscillator is 
tuned to the Larmor frequency, which results in 
depolarization of the atoms. The atoms attempt to equally 
repopulate the Zeeman energy levels. However, the helium 


lamp is still beaming polarized light energy into the 
absorbtion cell, causing the atoms to absorb light’ energy 
and rise to an excited energy levbel. This absorbtion of 
light energy is detected through the use of an infrared 
detector. The RF oscillator frequency producing maximum 
light absorbtion is called the resonant frequency, and is 
determined through the use of a servo loop from the infrared 
detector to the RF variable frequency oscillator. 
Therefore, any change in the Earth's magnetic field 
intensity will result in a change in the Larmor frequency of 
the helium atoms in the helium absorbtion cell. This new 
Larmor frequency will be detected by the ASQ-81 
magnetometer. Since the gyromagnetic ratio of helium is 
28.024 HZ per gamma, this detection of the resonant 
frequency provides a measurement of the Earth's magnetic 


field intensity at any given time. A change in the Earth's 


Li 





Magnetic field intensity could signal the presence of a 
submerged submarine. 

The output resonant frequency developed by the 
magnetometer is converted to a proportional output voltage 
which is filtered through the Magnetic Anomaly Detection 
(MAD) bandpass filters for environmental noise reduction 
and utilized to drive a chart recorder for observation by 


an operator. [Ref. 6] 


B. TRANSFER FUNCTIONS 

Transfer DineEuonce LOD them AN/ASO-S1 filters were 
obtained from the manufacturer of the AN/ASQ-81 detecting 
set, Texas Instruments of Dallas, Texas. These transfer 
functions are listed in Appendix A and are in the form of 


H(s), that is, the frequency domain, or S domain, where 


S = jw. The s-domain representation for transfer functions 
is routinely utilized to express output system 
characteristics for given system inputs. As the S domain 


representation is not utilzed further in this discussion, it 
will not be further explained. 

As the output signal of the ASQ-81 magnetometer is 
filtered through a fixed high-pass system, then through a 
selectable low pass system and a selectable high pass system 
(as shown in Figure 3.2 below), the transfer funtions are 


listed in this order. 
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SELECTABLE SELECTABLE 










HIGH PASS LOW PASS 


Further discussion will be made of the selectable filters 


later. 
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ieee! Cone rte k ING MODELLING OF SYSTEMS 


A. SEQUENCE REPRESENTATION OF TIME FUNCTIONS 
faobenal Kepresentation 

A signal can be defined as a function that conveys 
information, generally about the state or behavior of a 
physical system. Although signals can be represented in 
many ways, the information conveyed by the signal is 
contained in a pattern of variations of some form. Signals 
are represented mathematically as functions of one or. more 
independent variables, one of the most common of which is 
time. 

The independent Variable of the mathematical 
representation of a signal may be continuous or discrete. 
Continuous time signals are signals that are defined over 
continualy values of time and are therefore represented by 
continuous-variabled functions. Discrete time signals are 
defined at discrete time intervals and are therefore 
represented by functions whose independent variable(s) take 
on discrete values only. Discrete-time signals are 
represented as sequences of numbers. [Ref. 7] 

In addition to the fact that the independent 
variables can be either continuous or discrete, the signal 


amplitude can be either continuous or discrete. Digital 


Zo 





Signals are those for which both time and amplitude are 
discrete. Analog signals are those for which both time and 
amplitude are continuous. 

Digital signal processing deals with transformations 
of signals that are discrete in both time and amplitude, 
usually represented by sequences of numbers. The nth number 
in the sequence x being processed is usually represented as 
hem, and is formally written as: 

c— 1“ ie —CO< n < +o 
In general, an arbitrary sequence can be expressed as 
oe) 
nee x(k) d(n-k) 
k= -@& 
where d(n-k) is the unit sample at time k. In other words, 
an arbitrary sequence may be expressed as a sum of. scaled, 
shifted unit samples, where the scaling factor is equal to 
the amplitude of the sequence at that time. 
2. Linear Shift-Invariant Systems 

A system is defined mathematically as a unique 
transformation or operator that maps an input sequence 
[x(n)] into an output sequence [y(n)]. This is denoted 
as: 

y(n) = T[x(n) ] 


and is often depicted as in Figure 4.1. 
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Peeguire 4,1: Representation of Transformation of an Input 
Sequence to an Output Sequence. [Ref.7]} 
Classes of discrete time systems are defined by 
placing constraints on the transformation T{]}. 
The class of linear systems is defined by the 
Be@nciple of superposition. If y (n) and y (n) are the 
1 Z 
responses when x (n) and x (n) are the inputs, then a system 
] 
is linear if 
Rlaxe¢n) + "bx (n) |) = allx (mn) ] + bIix (n)] 
l 2 i ie 
wea en eeeee 1D nl) 
i Z 
for any arbitrary constants a and b. This, together with 
the concept of representing a sequence by a sum of delayed 
and scaled unit-sample sequences, suggests that a linear 
System can be characterized by its unit-sample response. 
Specifically, let h (n) be the response of the system _ to 


k 
d(n-k), a unit sample occurring at n=k. Then 


CO 


y COE | \ Ck eedCn= i) | or, 
La 


k= -©O 


ee 





y(n) = x(k) T[d(n-k)] = x(k ) oe 
k= -—OO k= -OO 
Thus the system response can be expressed in terms of the 
response of the system to d(n-k). 

The class of shift invariant systems is characterized 
by the property that if y(n) is the response to x(n), then 
y(n-k) is the response to x(n-k), where k is a positive or 
negative integer. When the index n is associated with time, 
shift-invariance corresponds’ to time-invariance. The 
property of shift invariance implies that if h(n) is the 
response to d(n), then the response to d(n-k) is’ simply 


h(n-k). Therefore BS 


y(n) = x(k) h(n-k) 
k= -—OO 
and any linear shift-invariant system is completely 
characterized by its unit-sample response h(n). 

A subclass of linear shift-invariant systems consists 
of those systems for which the input x(n) and the output 
y(n) satisfy an Nth-order linear constant-coefficient 
difference equation of the form 

N M 


> a y(n-k) = > b x(n-r) 
k r 


k=O r=0 





If the assumption of causality is made about’ the 
system, a linear difference equation provides an explicit 
relationship between the input to the system and the output 
of the system. This can be seen by rewriting the previous 


equation as 


N M 
y(n) = ». c y(n-k) + » d x(n-r) 
k ie 
k=] r=0 
where c =-a /a and c =b /a 
k le SO r rr 0 


Thus the nth value of the output can be computed from 
the nth value of the input and the N and M past values of 
the output and input, respectively. The difference equation 
not only represents the system for theoretical purposes, 
but it may also serve as a computational realization of the 
system. The z-Transform makes use of this property to 


realize systems. 


B. THE z-TRANSFORM 


il. Description of the z-Transform 





The z-transform plays an important role in the 


analysis and representation of discrete-time linear shift- 


invariant systems. The z-transform, X(z), of a sequence 
x(n) is defined as CO 
-n 
GZ es = Raa!) Zi 
n= -—- Cc 
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where z is a complex variable. This representation of the 
z-transform is referred to as the two-sided z transform. The 
one sided z-transform consists of the same summation for 
terms of n greater than or equal to Zero. For the case 
that x(n)=0 for n<O, the one sided and two sided 2z 
transforms are equivalent. 

By expressing the complex variable z in polar form as 

|W 
Z = a , the z-transform can be interpreted as the Fourier 
transform of x(n) multiplied by an exponential sequence. 
For r = 1, that is, for Izi = 1, the z-transform is equal to 
the Fourier tranform of the sequence. 

2. The Bilinear Transformation 
The transfer functions of analog systems are most 


often expressed in terms of s = jw (see section III B.). 


This corresponds to the analog frequency response of the 


system. This analog frequency response can be "mapped", 
that is, tranformed to the z-plane from the s-plane through 
the use of the bilinear transformation. The effect of 


utilizing the bilinear transformation is to convert a system 
transfer function in terms of the variable S into the system 
transfer function in terms of the variable Bo . The 


transformation itself is: 


and 


ILS) 





(2/T) + s 


where T is the sampling period, that is, the time between 
data samples. 

Thus a transform can be made from one plane to. the 
other. itn ocmmaniner,esthe transter function, H(z), of a 
system may be obtained, 

The bilinear tranformation equations may be shown to 
hold in general, and the use of this transformation may be 
shown to yield stable digital filters from stable analog 
filters [Ref. 7]. tine bilinear transformation maps the 
imaginary jw axis in the s-plane onto a unit circle (of the 
region of convergence) in the z-—-plane, with the left half 
s-plane mapped onto the region inside the circle and _ the 
right hand (region of instability) s-plane mapped onto’ the 
region outside this circle [Ref. 8]. A complete discussion 
of the z-transform is available in several texts, some of 


which are listed in the Bibliography. 


C. THE DIGITAL COMPUTATIONAL ALGORITHM 

In implementing a digital filter on a digital computer 
such as the IBM 3033, the input-output relationship of the 
signals through the system being synthesized must be 
converted to a computational algorithn. The algorithm is 
specified in terms of a set of basic computations of 


elements, For the implementation of discrete-time systems 
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described by linear constant Coerriclent difference 
equations, such as the AN/ASQ-81, it is convenient to choose 
as these elements the basic operations of addition, delay, 
and multiplication by a constant. The computational 
algorithm for implementing the filter is then defined by a 


structure or network consisting of an interconnection of 


these basic operations. For a system transfer function of 
the form M 
2" 
b z 
k 
k=0 Vez) 
H(z) = ----------- = —---- 
N Czy) 
-k 
l- » a z 
k 
eal 


the difference equation relating input and output is easily 
written down directly from the system function and is given 

by N M 
yCn) = | a y(n —~ k) + | b x(n - k) PReteen 7 | 

k k 
ee k=0 

This difference equation can be interpreted directly as a 
computational algorithm in which the delayed values of the 


input are multiplied by the coefficients b , the delayed 


k 
values of the output are multiplied by the coefficients 
a , and the resulting products are added. It is now easy 
k 
to see the process to be followed in obtaining the 


computational algorithm for the AN/ASQ-81 magnetometer 
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Eeanster function. The z-transform of the system transfer 
function is obtained through the use of the bilinear 
transformation, and is then converted into a difference 
equation relating input and output signals, thence to a 
FORTRAN computer program. A table of z-transforms of system 
functions is included in Appendix B. 

In the FORTRAN computer program realization of the total 
System computational algorithn, each filter block is 
transformed it © a separate difference equation and 
algorithn. This was done to enable a "building block" type 
approach to the program, and to minimize computational and 


moundoff errors. 


D. THE CASCADE FORM OF THE COMPUTATIONAL ALGORITHM 

Even though the direct form realization of the digital 
filter design may be perfectly satisfactory in a theoretical 
sense, it may be less than desirable in the context of 
realization through the use of a general purpose computer of 
fixed register length. The parameters of a digital filter 
are usually obtained with a high degree of accuracy, which 
results in a faithful realization of the desired system. 
When these parameters are quantized, as in a finite memory 
register within a computer, the frequency response of the 
resulting digital filter may differ appreciably from _ the 


original design. In fact, the quantized filter may fail to 
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meet design specifications although the unquantized filter 
does. [Ref. 7] 

imes, Sensitivity of the filter response to errors in the 
filter parameters is dependent upon the structure of the 
j[mereer realization. Therefore, in the event of an 
unacceptable change in the frequency response of the filter 
due to quantization errors, it is often possible to minimize 
the effect of these errors through an alternate filter 
realization structure. An alternate structure to the 
previoulsy discussed direct form realization is the cascade 
morm realization. 

The direct form network structures were obtained directly 
from the system function H(z) written in the form of a ratio 
of sums. iitas ratlourts factored into a4) product of 
polynomials of the form 

[(N4+1)/2] _} b 


H(z) = A | [| aaa a--------------------- 


a Z Z 
kel 1k 2k 


this product represents a general distribution of poles’ and 
zeros and suggests ae set of structures consisting of a 
cascade of first and second-order subsystems. There is 
considerable freedom in the design of the subsystems, but it 
is best to realize the systems using a minimum of storage. 
The expression of H(z) in this form indicates’ the 


presence of poles and zeros in pairs. If poles and zeros 
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are not present in pairs, one of the coefficients B or a 
aks 2k 
will be zero as appropriate. An implementation of such a 
cascade structure with the use of minimum memory can _ be 
obtained through a direct form II realization of each second 
order subsystem using techniques similar to the direct form 
implementation utilized previously. A cascade realization 


of a sixth-order system, such as the ASQ-81 system, using a 


direct form II realization of each second order’ subsystem 


would appear as in Figure 4.2 below. 





Figure 4.2: Cascade Structure With a Direct Form II 
Realization of Each Second Order Subsystem. 

Ref. 7] 
There is, theoretically, considerable flexibility in the 


manner in which the poles and zeros are paired together and 
in the order in which the resulting second-order subsystems 
are cascaded. However, although alli such pairings and 
orderings are equivalent for infinite-precision arithmetic, 
they may differ considerably in practice owing to finite 


word length effects of roundoff and truncation. 
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Figure 5.1: Output of First Stage Filter of Digital Filter 
Computer Program With Sinusoidal Input in Simulation. 
Unfortunately, the second stage output of the filter 

showed an instability within the program design, indicated 

by the output of the filter being a sinusoid of increasing 


magnitude, as indicated in Figure 5.2 below. 
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Figure 5.2: Output of Second Stage Filter Design With Input 
Ofrm@a Sinweeud . 
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The stability of the third stage of the filter design 
was investigated by inputting the sinusoid directly to the 
Meaird filter, and found to be stable. A check of the 
derivation of the equations, coefficients, and programming 
steps of the second (unstable) filter of the design failed 
to indicate the cause of the instability. 

Computation of the poles of the z transfer function, 
H(z), of the second stage of the filter confirmed the 
instability of the design. The poles were computed to be: 
Seeoet O-1218i, 1.07 + 0.13401, 0.8611, and 1.1564. Of 
these six poles, three lie outside the region of convergence 
for the z-plane, that is, within the unit circle discussed 
previously in Chapter IV. 

The second stage of the filter was therefore 
redesigned using the cascade form of the direct’ form 
realization (direct form II), and tested in simulation. A 
copy of the software used in the simulation is enclosed in 
Appendix F. 

The output of all three filter stages of the program 
were stable, as indicated in Figures 5.3 through 5.7 below. 
The amplitude decrease and phase shift expected were 
observed. The "damped overshoot" of the second stage output 
is due to the fact that, for values of the input function 
prior to time zero in the simulation, utilized in the input- 
output signal difference equations for the filter, the input 


Signal was set at Q. This resulted in an instantaneous 
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change of the input signal from 0 to the finite value 
introduced in the simulation at time 0+. The “overshoot” of 
Phe: filter is the filter's attempt to "match" this 
instantaneous jump in magnitude of the input signal. When 
the input signal to the filter in the simulation is Zero 


at time zero, this overshoot effect does not occur. 
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mgure 3.3: Input Signal to Digital Filter Program. A 
Sinusoid of Frequency 0.1 HZ and Amplitude + 1. 
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Pagure 5.4: Output of First Stage of Digital Filter Program. 
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Figure 5.5: Output of Second Stage of Digital Filter 
Program. 
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Meemre 5.6: Output of Third Stage of Digital Filter Progran. 


This simulation was run with inputs of sinusoids’ of 
various frequencies in order to check the stability of the 
filter design at frequencies throughout the operating range 
of the AN/ASQ-81 magnetometer. In all cases, the design was 
stable, and the expected amplitude changes and phase shifts 
occured. 

2. Noiselike Inputs 

The simulation was also run with inputs consisting of 
a sinusoid of a frequency which should be passed through the 
AN/ASQ-81 added to sinusoids of frequencies which should 


have been filtered by the magnetometer and random noise. The 


filter performed as expected, with the sinusoid of a 
passable frequency passed by the filter, and spurious noise 
and sinusoids attenuated severely. The results of a 


Simulation consisting of a sinusoid of passable frequency, a 
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filterable sinusoid, and uniformly distributed random noise, 


all of amplitude + 1, are presented in Figures 5.7/7 through 


5.10 below. 
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Meemre o./: Input to Filter - 0.1 HZ Sinusoid, 10 HZ 
Sinusoid, Uniformly Distributed Random Noise of Amplitude 
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FIctwemseos Output of First Filter Stage. 
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Rigunwe S210: Final Filter Stage Output. 
As can be seen, the digital filter program succeeds 


in filtering out random noise and signals of frequency 
components above the band pass of the magnetometer. 
In order to ensure that the digital filter 


representation of the magnetometer has the same amplitude 
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versus frequency Ghawacteristi¢es of the AN/ASQ-81 
magnetometer, a simulation program was written which inputs 
sinusoids of varying frequencies and computes the Root Mean 
Square (RMS) value of the filter output and the signal 
input, then computes the decibel (dB) attenuation of the 
filter at that frequency. AEECODY. OL thas program. is 
included in Appendix G. A plot was made of the dB 


attenuation versus frequency for the filter and compared 
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meeoure 5,11: Plot of Attenuation Versus Frequency for 


Sinusoidal Inputs for Digital Filter and Analog Filter. 
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with the measured frequency performance of the AN/ASQ-81 
magnetometer, which was supplied by Texas Instruments, Inc., 
and does not include the effects of the fixed high pass 
Pieter , Consequently, the data shown in Figure 5.1ll is a 
comparison of the data supplied by Texas Instruments and the 
output of the test program, which also does not include the 
fixed high pass filter. As can be seen, the performance of 
the filter is extremely similar to that of the magnetometer 


itself. 


3. Anderson Function Simulations 

The next step in the simulation phase was’ the 
introduction to the filter of Anderson function simulations. 
The shape of the signal amplitude of the output of a 
Magnetometer passing through the sphere of influence of a 
magnetic anomaly (submarine) is a function of the dip angle 
of the geomagnetic field, the magnetic heading of the track 
of the magnetometer (or the aircraft), the Magnetic 
heading of the anomaly (submarine) dipole, and the 
lateral range between the magnetometer (aircraft) and the 
anomaly. Anderson functions [Ref. 9] are mathematical 
representations of three basic components of signals’ which, 
when taken in various linear combinations, describe the 
shape of these anomaiy signals. The equations for’ the 


Anderson functions are: 
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3 3s) 


/eirst Anderson Function) f Glasby.) 


(velocity)x(time) 


range at CPA 


where B 


or, a dimensionless parameter. defined as the distance 
traveled along the magnetometer (aircraft) track divided by 
the slant range at closestpoint of approach (CPA) 
(Second Anderson Function) f= Bx £ 
(Third Anderson Function) _ = Bx — = * x f 
2 l 0 

The Anderson functions were introduced into the 
filter program in a noise-free signal environment in order 
to observe the output signal and ensure that it was a "MAD- 
like" signal. A rigorous determination of the actual output 


Signal would have been extremely difficult to obtain, so a 


comparison was made with the output of a computer simulation 


program Drovaded to SNPs. (by) Mr. Jige- Rice ot Texas 
Instruments. When the sampling rate of the program was 
adjusted to equal that of the Texas Instruments program, 8 
HZ, the two program outputs were observed to be very 
Samilar. The Anderson function simulation inputs.) and 


outputs of the program are depicted in Figures 5.12 through 
el. The Texas Instruments program outputs were obtained 
in the form of time series plots of discontinuous’ data 
points, and were therefore not conducive to replotting for 


comparison. 
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meewre 3,12: First Anderson Function Input. CPA at Time 
2.5 Minutes. 
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puewre 5.13: Filter Output for First Anderson Function 
np ut 
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Figure 5.15; Filter Output for Second Anderson Function. 
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Figure 5.16: Third Anderson Function Input. CPA at Time 
2.5 Minutes. 
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Heiter .t/seraiiter Output tor Third Anderson Function. 


The filter output for all three Anderson function 
inputs did appear to be "MAD-like" signals, and did closely 
resemble the Simulation output obtained from Texas 
Instruments, Inc. 

4. Impulse Function Response 
The response of the filter program was also observed 


When the input was a unit impulse function. Again, the 
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output was compared to that of the Texas Instruments' 
computer program. The outputs of the two programs were 
observed to be, again, very similar, as can be seen in 
Figure 5.18, where the response of the NPGS filter is 
represented by a solid line and that of the Texas 


Instruments filter by a chain-dash line. The abrupt "jumps" 
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of the Texas Instruments response are due to the translation 
of the output plot supplied to this plot. The plot supplied 
by Texas Instruments was, again, discontinuous points of 


poor resolution, and it was necessary to interpolate values 
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in order to generate Figure 5.18. fis resulted in the 
broken appearance of the plot. DVetmcommcnoms iii Laricy of 


the outputs can be observed. 


B. EQUIPMENT SETUP 

Following the simulation phase of the experiment, actual 
magnetic field measurements were introduced to the filter in 
order to test the response of the filter. Magnetic field 
measurements were made at the La Mesa field test site near 
the Naval Postgraduate School in Monterey. The output of an 
AN/ASQ-81 magnetometer, a Schonstedt magnetic field sensor, 
and the school's coil sensor, oriented along the Earth's 
magnetic field, were pulse code modulated (PCM) and 
transmitted via VHF radio to recording devices at the 
Postgraduate school. The recording of a two hour long data 
collection period was transferred to digital data tape for 
use by the school's IBM3033 general purpose mainframe 
computer. 

In the first test of the digital filter program, the 
output of the Schonstedt sensor, which represents 
fluctuations of the Earth's total field, was used as'- the 
input to the computer program. A comparison of the output 
of the computer program, with this approximation to’ the 
total field fluctuations as input, to the output of the 
AN/ASQ-81 should provide an indication of the’ proper 


Pimcrioning Of the computer filter progran. The results of 
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the test are shown in Figures 5.19 through 5.21 on the pages 
following. Figure 5.19, the Schonstedt sensor output, shows 
several instances of PCM dropouts, that is, occasions where 
the pulse code modulation signal was not correctly read by 
the computer for some reason. At such occurances, the data 
point value used by the computer is a random number and does 
not reflect the true value of the data. The problem with 
these PCM dropouts is that the computer does not recognize 
them as invalid data points and will use them 10 
computations. This can (and does) cause problems in the 
computation of Fourier transforms, spectral characteristics, 
etc. Additionally, this will also impact the proper 
functioning of the digital filter program which is’ the 
subject of this thesis. PCM dropouts are visible at times 
6, 8, 142, and 220 through 226 seconds on the plot of the 
schonstedt sensor output. An examination of Figure 5.20, 
the filter program output, reveals the programs attempt to 
"follow" these PCM dropouts. It should be recalled that 
previous simulations indicated the filter's tendency to 
"follow" sudden changes in the input’ signal, with a 
relaxation time required for the filter to steady out. This 
effect is apparent in the output of the filter program at 
times corresponding to those of the PCM dropouts in the 
Schonstedt sensor's time series plot. It can be seen that 
this overshoot tendency resulted in an output significantly 


different from the actual AN/ASQ-81 output at these times. 
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If the PCM dropout induced differences are neglected, it can 
be seen that the shape of the output of the filter program 
is remarkably similar to that of the AN/ASQ-81 magnetometer, 
although noisier. Note that the output of the AN/ASQ-81 
magnetometer exceded the maximum voltage amplitude which the 
pre-amplifiers of the data collection system were able to 
handle and resulted in a truncated signal from time 40 to 
time 60 seconds. It can still be seen, however, that the 
filter program output is very similar to the time signal 
which would have been displayed without this truncation. 

It should be noted that the amplitudes of the time series 
Signals of the program output and the AN/ASQ-81 magnetometer 
differ considerably. In the case of the program output, the 
peak amplitudes are on the order of 1.4 nanoteslas along the 
vertical scale, while the peak amplitudes of the output of 
the AN/ASQ-81 magnetometer are on the order of Oey 
nanoteslas along the vertical scale. This is because’ the 
input signal to the filter program is an approximation to 
the total field difference time series signal, and _ some 
amplitude difference could reasonably be expected. The 
intent of this initial test was to investigate the output 
time series’7 shape, and an exact correlation was not 
expected. It is worth noting that the digital filter 
program will perform its function on any time series signal, 


regardless of units. This means that a signal may be 


Sh 





operated upon either before or after conversion from 
whatever units it was originally measured to magnetic field 
strength units. 

Therefore it appears that the digital filter program is 
functioning properly. When a close approximation to the 
fluctuations of the total field time series signal is used 
as the input to the computer program, the output of the 
program is similar to the time series output of an AN/ASQ-81 
Magnetometer. 

The final stage in the testing process was a conversion 
of the time series output voltage signal of the coil antenna 
sensor, which was aligned along the Earth's magnetic field, 
into a total field fluctuation time series representation 
for the same time period as before, and then to use this as 
the input to the digital filter program. A comparison of 
the resultant time series output of the program with the 
actual AN/ASQ-81 magnetometer output would validate the 
proper functioning of the program. 

Conversion of the time series antenna sensor output 
voltage signal into total field fluctuations in nanoteslas 
was accomplished through the use of a computer program 
designed by Capt. Kurt Stevens, USAF, a student at the Naval 
Postgraduate School, as his Master's thesis [Ref 10]. The 
output voltage time series is stored in an array, then a 
Fourier transform is performed on the stored data, resulting 


in the Fourier spectrum of the data. This spectrum is 
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corrected for the characteristics of the coil antenna sensor 
to obtain the Fourier spectrum of the total field data. A 
reverse Fourier transform gives the time series signal for 
total magnetic field in nanoteslas. 

This time series signal was used as the input to the 
digital filter program and compared with the output of the 
AN/ASQ-81 magnetometer. Figures 5.22 through 5.24 show the 
raw coil antenna data, the total field time series data, and 
the program output time series for a 6 minute period of the 
test. Figure 5.22 shows the raw coil antenna data series. 
The number of PCM dropouts should be noted, as these will 
influence the performance of the filter program. Figure 
5.23 shows the computed total field time series. Note that 
the PCM dropouts evident on the raw time series plot are 
evident on the computed total field time series plot also, 
and thus inputted to the filter program as valid data 
points. Additionally, there are two "jumps" in the plot of 
total field fluctuation (Figure 5.23) which are also 
inputted to the filter program as valid data points. These 
"jumps" are located at 128 and 256 seconds and are caused by 
the method of processing blocks of data for the conversion 
Bo total field fluctuation, A block of 128 seconds of data 
is processed at a time, and the results of each block are 
Stored in an array. This results in a slight amplitude 


difference between the last data point of one block and _ the 
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Mmimest data point of the next block of data. ihre — slight 
difference is manifested as a signal jump. 

A comparison of Figures 5.24 (program output) and 5.21 
(AN/ASQ-81 output) show that the filter program gives a time 
series output very similar to that of the actual 
magnetometer. The first 20 seconds of the program output is 
somewhat dissimilar to that of the AN/ASQ-81, due either to 
the initial "start up" delay of the filter program or _ to 
distortion of the total field fluctuation time series. 
There is a PCM dropout at time 11 seconds which contributed 
moethe distortion. 

Following this, however, it can be seen that the program 
output is very similar to that of the magnetometer, except 
at 128 and 256 seconds, which show the effects of the false 
signal jumps caused by the total field PwC E Wat Von 
conversion. There is also a noticeable four to five second 
time delay between the AN/ASQ-81 output and that of the 
miter program. As this delay is not evident in a 
comparison of the AN/ASQ-81 output and that of the filter 
program with the Schonstedt sensor as the input, it can be 
inferred that this time delay is caused either by the 
program which converts the raw coil data to total field 
fluctuation data, or by a phase (and hence time) change of 
the voltage signal due to the coil sensor itself. A 
comparison of the raw coil data in Figure 5.22 to the 


Souverted coil data in Figure 5.23 indicates no time shift, 


a7 





and hence the deduction can be made that there is a time 
delay inherent within the coil sensor itself. 

Other than the differences of the four to five second 
time delay and the distortions caused by the false signal 
jumps, the program output is extremely similar to the output 


of the AN/ASQ-81 magnetometer. 
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VI. CONCLUSIONS 

The intent of this thesis was to design and test a 
digital filter computer program which would, when given a 
time series input of fluctuations in the total magnetic 
field, deliver an output time series representation of the 
output of an AN/ASQ-81 magnetometer. This purpose has been 
realized. 

The computer program contained in Appendix [I has. been 
proven to output a time series signal which is very similar 
to that of the magnetometer. The major limitations of the 
output signal are a finite time delay of about five seconds 
between the AN/ASQ-81 magnetometer signal and the output 
signal of the program, a sensitivity of the program to false 
data points such as those caused by PCM dropouts and _ false 
signal jumps caused by processing large data blocks, and the 
inherent limitations of the program caused by its dependence 
on the use of digital data tapes and the IBM 3033 mainframe 
computer. 

The five second time delay is not considered to be an 
important limitation to the program, as it was intended as a 
research tool for programs currently in progress at the 
Naval Postgraduate School. Instances where this time delay 
might become important would be in areas of simultaneous 
comparison of the program output signal with an actual 


mMagnetometer, in target location algorithms using time 


De, 





delays, or in correlation studies between different sensors. 
fae correlation studies using coil sensors, the effects of 
the time delay would cancel out, as all coil outputs would 
be similarly delayed. In target location algorithms, the 
target location errors due to the time delay could be 
adjusted for simply while in computer simulation, and flight 
testing could not reasonably be accomplished without the use 
of an actual magnetometer as the sensor. Lastly, in a 
comparison of the program output with an actual sensor, the 
time delay can, again, be compensated for. In short, these 
limitations are not considered excessive, especially as the 
apparent cause for the delay is not the filter program. 

In the primarily intended purpose of the filter progran, 
magnetic noise studies, the time delay is not considered to 
be a problen. 

The problem of false data points caused by PCM dropouts 
and Signal jumps due to conversion. to total field 
fluctuations is more _ serious. False data points cause 
inaccuracies in the output time series and could adversely 
affect later projects. Unfortunately the PCM dropout 
problem is one which is endemic to the data collection 
system presently being used at the postgraduate school, and 
not to the filter program itself. It is imperative that 
users of this program are aware of the PCM dropout problem 


and of the effects it may entail upon their specific 





research. A large number of PCM dropouts in atime series 
could render that series unusable. Similarly the case of 
the false data jumps caused by conversion to totai field 
fluctuations is not within the filter program. Poet ne 
investigation of this problem is necessary in order to 
Siminate it. 

The last problem, that of reliance upon the digital data 
tape/IBM 3033 computer system, is, like the PCM problem, one 
which is not endimic to the filter program but rather to the 
data collection system being used. A change of data 
collection system may, at some future time, remove’ the 
reliance upon the PCM/digital tape/IBM 3033 data system (and 
hence too the data block conversion problem which results in 
falso data jumps), but this is unlikely at this time. Users 
should be aware of this dependence and of possible effects 


upon specific research projects. 
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APPENDIX A 


AN/ASQ-81 FILTER TRANSFER FUNCTIONS 


Fixed High Pass Transfer Function: 


80 S- + 20 S$ + 1 


Selectable High Pass Transfer Functions: 


A. 0.04 HZ H(S) = 
2 2 
40.82834 § 45.28317 § 
----------- GR One SST SSS SSS S= x Son Sg SH SS SSS SS SSS SS 
mmeoe34 $- + 12.52096 S +1 45.28317 S~ + 11.00999 § +1 
2 

So) SUBIR) 

x< == om oo of © © © © © © © © o © © © © © we om © oe we ww ww ww = w= = 
Bs 6688 S- + 7.41498 S + 1 
B 0.06 HZ WS) < 
2 2 
Tee aseues 20.12587 § 
oo om we om oe © oo ow ow wo ow = om ow om ow = = oo & Oo = = = ve San to Ss Sa oS = =e 
Memiaso1 S° «+ 8.34727 § + 1 20.12587 S- + 7.33999 $ +1 
2 

25.58964 § 

Xo ----------5----- -- --- -- -- -- -- 
PEMSCO64 Gree 4.94332 5 < 1 
oe 0.08 HZ H(S) = 
2 ” 
10.20708 § 11.32080 § 
MOM20708 S$- + 6.26045 $ + 1 11.32080 S° + 5.50500 §$ + 1 
2 

14.39417 § 

xX me ge = © om Oo Oo o® Om w= = oO Oo = Of = Oo OF = OO OO Ow ws = w= = 


WOR 7 S- + 3.707495 + 1 
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2 
6.53253 § 7 DWORSA S 
werner nr -— imine ee x wis 1S in ee oS 
feo253)«6|CUSlC(‘it+#«CO55«wL 00836 S + 1 7.24531 S + 4.40400 S + 1 
2 
5) DIPD| S 
x Cee eee eS oo Se SS sa OS 9 nn nn aa a 
9.21227 $* + 2.96599 S$ +1 

Selectable Low Pass Transfer Functions 
fee 0.2 HZ HCS = 

1 1 
me ee ee ew es we ws se es we ww ee ee ee = = x Rn 
fees S24 1.0741 S$ + 1 COS OlmsS= 6.6200 S 4 1 
B. 0.4 HZ Hasn = 

1 1 
|= © ow © © © ow © © © © Ow © OO © © 8 = Ow = ow = @® oo = oe De Sak. Get AG a eee seine 
0.07858 S* + 0.53706 $ + 1 OnOGI5 2S ao 31044 9 4 I 
0.6 HZ CS) a 

1 1 
as?) ~~~ ---------—— X --------- We = Seat ee 
mmee492 S~ + 0.35804 $ + 1 GO) O2mmemc = Ouaeiccbe.S + 1 
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Me Je jet 1D Ip 5) 


AN/ASQ-81 Z TRANSFORM FILTER TRANSFER FUNCTIONS FOR DIRECT 
FORM I REALIZATION 


For fixed high pass filter: 





1 - AFHP1*Z-” - AFHP2*Z 


where BFHPO, BFHP1, BFHP2, AFHP1, . AFHP2 are constants 
tabulated in Appendix D. 


For selectable high pass filter: 
— = =) -3 
BSHPO + BSHP1*Z + BSHP2*Z + BSHP3*Z 


i 


+BSHP4*Z~* + BSHP5*Z > + BSHP6*Z ~° 


AP ASHP1*Z 1 = ASHP2*Z ~- ~ ASHP3*Z 


4 5 


- ASHP4*Z ° - ASHP5*Z ~ - ASHP6*Z~ 


where, for low frequency cutoff of 0.04 HZ: 


BSHPO = 0799471378 

BoHPl = -5.9682827 ASHP] = 5.9894021 
Bemte2 = 14.920707 ASHP2 = -14.94/7051 
BSHP3 = -19.894276 ASHP3 = 19.894225 
BSHP4 = 14.920707 ASHP4S = -14.894327 
BSHP5 = -5.9682817 ASHP5 = 3.9472141 
BSHP6 = 0.99471372 ASHP6 = -0.98945296 
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For low frequency cutoff of 0.06 HZ: 


BSHPO 
BoP | 
BeHe 2 
PouP 3 
BSHP4 
BSHP5 


BSHP6 


For low 


BSHPO 
Bee 1 
Bo P 2 
BoP 3 
BSHP4 
BSHPS 


BoHP 6 


For low 


BSHPO 
Beae 1 
BoP 2 
BoHP 3 
BSHP4 
BSHP5 


BSHP6 


On 99 20E 15 


=5 9524926 


14 
-19 
14 


=9,9524927 


OE 20821)2 


1681252 
~841643 


mo oe 


ASHP! 
ASHP2 
Noes Beas) 
ASHP4 
ASHP 5 


ASHP6 


frequency cutoff of 0.08 HZ; 


0.98945806 


-5 
14 
-19 
14 


-5 


0.98945802 


Frequency cutoff of O.1 
0.98684156 


-5.9210493 


-9367483 
.841871 
WSSU Il 
~841871 


~9367476 


14.,802623 


-19 
14 


-5 


0.9868415 


272/06 31 
~o02623 


~9210491 
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ASHP 1 
ASHP 2 
ASHP 3 
ASHP4 
ASHP 5 


ASHP 6 


Lae 


ASHP 1 
ASHP 2 
ASHP3 
ASHP4 
ASHP = 


ASHP6 


-14, 


OS, 


-14, 


-14, 
gh 
-14, 


-14, 
19. 


EAP 


-9841070 


920650 


1241528 


841757 


SPOS SIS, 
~98422128 


.9788144 


894276 
UNS) SS) 
feos 


~-8948841 


eo 2027 Z0 


~9735244 


867941 
7/015 16 


737149 


. 8688889 


~97384065 





For selectable low pass filter: 
-l -2 -3 me 
BSLPO + BSLP1*Z + BSLP2*Z + BSLP3*Z + BSLP4*Z 





1 - ASLP1*Z71 — ASLP2*Z ~* - ASLP3*Z7> — ASLP4*z ~* 
where BSLPO, BSLP1, BSLP2, BSLP3, BSLP4, ASLP1, ASLP2, 


ASLP3, ASLP4 are constants tabulated in Appendix D. 
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APPENDIX C 


AN/ASQ-81 2Z TRANSFORM FILTER TRANSFER FUNCTIONS DIRECT 
FORM II REALIZATION 


For fixed high pass filter: 





1 - AFHP1*zZ~+ - AFHP2*Z~? 
where BFHPO, BFHP1, BFHP2, AFHPI1, AFHP2 are constants 


tabulated in Appendix D. 


For selectable high pass filter: 





-1 —2 
1 — 2*Z + Z 
fee) = ASHPI] x eee ge eee ge X 
1 - ASHP3*Z - ASHP2*Z 
i =) -1 -2 
Pe— 2*7, + Z 1 - 2*Z + Z 
—.. Se ee °s 2) 2 eee 9 ee D 
1 - ASHP4*Z - ASHP5*Z P= BAS HE Oe, -~ ASHP7*Z 


mere ASHP!, ASHP2, ASHP3, ASHP4, ASHP5, ASHP6, ASHP? are 
constants and tabulated in Appendix D. 
For selectable low pass filter: 


=I = =e =) 
BSLPO + BSLP1*Z + BSLP2*Z + BSLP3*Z + BSLP4*Z 





1 — ASLP1*Z > - ASLP2¥*Z~° - ASLP3*Z °- ASLP4*z* 


PeemesboLrPO, BSLPi, BSLP2, BSLP3, BSLP4, ASLP1, ASLP2, 


ASLP3, ASLP4 are constants tabulated in Appendix D. 
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APPENDIX D 


Z TRANSFORM REALIZATION DIFFERENCE EQUATIONS 


With reference to Figures 3.2 and 4.2, the following 
difference equations are used to model the AN/ASQ-81 
magnetometer filter transfer functions. The input to the 
fixed high pass filter is called SIG(I), where I is’ the 
current data sample. The output of the fixed high pass 
filter, which is the input to the selectable high pass 
filter, is YO(I), and the output of the selectable high 
pass filter, the input to the selectable low pass filter, 
meee called YPOCI). The output of the filter is called 
ASQ(I). (I-1) denotes a time delay of one sample, and so 
forth, and the symbol * denotes multiplication. 

For the fixed high pass filter: 


YO(I1)=BFHPO*SIG(I) + BFHP1*SIG(I-1) + BFHP2*SIG(I-2) + 
AFHP1*YO(I-1) + AFHP2*YO(I-2) 


where: 

BFHPO = 0.9980499222938581 

BFHP1 = -1.9960998445877161 AB Sl IIOOF Sa 7 16543922 
BFHP2 = 0.998049922238581 AFHP2 = -0.99610136/4910398 


For the selectable high pass filters: 
XI(I) = ASHP1*YO(I) + ASHP2*XI(I-2) + ASHP3*XI(I-3) 
XII¢(I) = XICI) + XI(I-2) - 2*xXI(CI-1) 
XITICI) = XIICI) + ASHP4*XIITICI-1) + ASHPS*XITI(CI-2) 


XIV(I) = XIIICI) - 2*XIITICI-1) + xXIII¢CI-2) 
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XV(I) = 


YPO(I) 


For the low frequency cutoff at 0.04 HZ: 


eoae | = 
ASHP3 = 
ASHP5 = 
mode 7 = 
For the 
pone] = 
ASHP3 = 
ASHP5 = 
ome 7 = 
For the 
aone l= 
ASHP3 = 
ASHP5 = 
ASHP7 = 
For the 
ASHPI = 
ASHP3 = 
ASHP5 = 
ASHP? = 
For the 


ASQ(T) 


selectable low pass filters: 


XIV(I) + ASHP6*XV(I-1) + ASHP7¥*XV(I-2) 


XV(I) 


0.994713789347288 
ioe eZ b3s/ 256322473 
-~0.9962082013015601 


-0.9979897681491607 


low frequency cutoff 


One 22087127 7199399 
ie 99262472455 17998 
-0.9943177045263504 


=0:.9909861717656565 


low frequency cutoff 


Or ISIS OUD D6G75767 
1.9904384565912725 
-0.9924307799269113 


-0.9959835860201267 


low frequency cutoff 


0.9869415560096681 
1.9880549117849344 
-0.9905474442556225 


-0.9949820174650785 


mee rycen) + XV (I=2) 


ASHP2 
ASHP4 


ASHP6 


ASHP 2 
ASHP4 


ASHP6 


=(). 
i. 


Ue 


MZ 


=) 


1 


1 


ateoncs HZ: 
-O. 


Le 


ASLP1*ASQ(I-1) + ASLP2*ASQ(I-2) 


ASLP4*ASQ(I-4) + BSLPO*YPO(I) + 


BSLP2*YPO(I-2) + BSLP3*YPO(I-3) 
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oy ZOO 01 57 S62 
9962028201103847 


9979855321466768 


99283813061 74365 


.9943056083414792 


-9969766455640039 


9904602761133 7722 


9924092959984783 


SSS Be OSs ers Ilsa 


BySoe 72701077 949) 
p70 5159074397893 
~9949555824611562 


+ ASLP3*ASQ(I-3) 
BSLP1*YPO(I-1) 


+ BSLP4*YPO(I-4) 





Bemeche high frequency cutoff at 0.2 HZ: 

BSLPO = 0.0000000452616229 

BSLP1 = 0.0000001810464917 AS ri = oe oUee aes 59 L077 
BSLP2 = 0.0000002715697375 ASLP2 = -5.7285022156249328 


BSLP3 = 0.0000001810464917 p Wey iy ae, ao homes LOUIS 


Boer4 = 0,.0000000452616229 ASLP4 = -0.9119356638511430 
For the high frequency cutoff at 0.4 HZ: 


BSLPO = 0.0000006918001209 


BSLP1 = 0.0000027672004837 ASLP1 8.01 73771378993420 


BSLP2 = 0.0000041508007256 ASLP2 -5.4670046844062743 


BSLP3 = 0.0000027672004837 ASLP3 = 3.4812554457127576 
BSLP4 = 0.0000006918001209 ASLP4 = -0.8316389680077603 
For the high frequency cutoff at 0.6 HZ: 

BSLPO = 0.0000033463317975 


BSLP1 = 0.0000133853271900 ASLP 1 3.7274299052305002 


BSLP2 = 0.0000200779907850 ASLP2 —5.2152772583906819 


EobPo = 0.0000133853271900 SE ge: 3.2462216641520829 


BSLP4 = 0.0000033463317975 ASLP4 —-0.7584278523006615 
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DIGITAL SOFTWARE FOR SIMULATION - DIRECT FORM (SINUSOIDS AS INPUT) 
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